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Abstract 

The aim of this discussion is to expose incorrect results in a previous IJSS article. 
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Preliminaries. The article by iRocha de Faria et al.l (120071 ) under discus 



sion is concerned with the evaluation of the perturbation undergone by the po- 
tential energy of a domain Q (in a 2-D, scalar Laplace equation setting) when 
a disk B £ of small radius e centered at a given location x G Q is removed from 
Q, assuming either Neumann or Dirichlet conditions on the boundary of the 
small 'hole' thus created. In each case, the potential energy ip(Q E ) of the punc- 
tured domain Q £ = Q \ B £ is expanded about e = so that the first two terms 
of the perturbation are given. The first (leading) term is the well-documented 
topological derivative of ip. The article under discussion places, logically, its 
main focus on the next term of the expansion. However, it contains incorrrect 
results, as shown in this discussion. In what follows, equations referenced with 
Arabic numbers refer to those of the article under discussion. 



Topological expa nsion: Neumann conditi on on the hole. In the main 
result proposed by IRocha de Faria et al.l (120071 ) for this case, namely expres- 
sion (37) for the topological expansion of the potential energy, the first term 
(whose order is 0(e 2 )) is correct but the second (whose order is 0(e 4 )) is not 
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as it lacks a contrib ution of the same order related to the external boundary 
(see iBonnetl . l2006al . for a similar study in 3-D linear acoustics). 



This error can be explained as follows. Equation (37) is based on an expansion 
of 

-1^(0.) = -± / (Vu £ .e e ) 2 ds (i) 

de 2 JdB e 

up to order 0(e 3 ) (where (e n e#) are the unit vectors associated with polar 
coordinates (r, 9) originating at the center of B £ ). Since ds = ed9 on dB £ , this 
task requires expanding { ^u F {x).eo) 2 to order 0(e 2 ) f or x 6 dB £ . The latter 



operation is carried out in iRocha de Faria et al.l (120071 ) by evaluating Vu e (a;) 
from the 0(e 2 ) expansion (23) of u £ . However, expansion (23) evaluated on 
dB £ gives 

Vu £ (x).e e = 2Vu(x).e e + 2eWu(x) : (e r <g> e e ) + 0(e 2 ) (x E dB £ ), 

and is therefore not suitable for expanding {Wu £ .ee) 2 to order 0(e 2 ) as it lacks 
the necessary 0(e 2 ) contribution to Vu £ .ee. The missing 0(e 2 ) term stems 
from the 0(e 3 ) contribution to u £ and is in fact non-local as it is expressed in 
terms of quantities on dfl rather than higher-order gradients of u at x. 

The incorrectness of result (37) can be further demonstrated on a simple 
analytical example. Consider the 2-D domain Q £ enclosed by two concentric 
circles of radii e and a, i.e. dB £ = {(r,8) r = e} and dQ — {(r, 9) r = a} in 
terms of polar coordinates (r, 9). The solution u £ of the Laplace equation with 
boundary conditions 

u, n = (r = e), u tn = q = cos 9 (r = a) 

and the corresponding reference solution u when there is no hole are respec- 
tively given (up to an arbitrary additive constant) by 

a 2 ( e 2 \ 

u e{ r -,9) — -t, ^ [r -\ ) cos6>, u(r,9)=rcos9 (ii) 

a 2 — e 2 \ r J 

Note that the reference solution u is such that Vn(i) = cos#e r — sin^eg and 
VVh(x) = 0. Then, a simple calculation gives 



ip(tt £ ) = - ( Vu £ .Vu £ dV - [ qu £ ds = ~ [ qu £ di 
2 Jn E Jan 2 Jan 

Expanding ip(Q £ ) to order 0(e 4 ) gives 



7ra 2 a 2 + e 2 
2 a 2 -e 2 



i J (n £ ) = -^--ne 2 --e i + o(e i ) (hi) 
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while equation (37) incorrectly gives the expansion as 



V>(ft 6 ) = -— - 7T6 2 - X 6 4 + 0(6 4 ) (iv) 

Note that the error in (iv) vanishes as dQ is rejected to infinity, i.e. as the 
influence of the external boundary goes away. This is analogous to secondary 
reflection effects in small-obstacle approximations for wave problems. 



Topological expansion: Dirichlet condition on the hole. The topo- 
logical expansion (38) is also not correct. Expansion (38) states that 



-i 

Loge 



[u(x)Y + 7r||Vu(iB)||V + o(e 2 ). 



However, another simple analytical example again allows to show that the 
second term in (v), is not correct. With the domain fl e defined as before, the 
solution u £ of the Laplace equation with boundary conditions 

u = (r = e), u — A (r — a) 

and the corresponding reference solution u are respectively given by 

Log (r/e) 



u e (r,6) = A- 

Log (o/e) 

The potential energy is therefore 



u(r,0) = A 



me 



2 Jn e 



VUr.VUr dV = ~ 



du £ 
2 Jan on 



7TA 2 



7TA 2 



Log(a/e) Log a — Log £ 



Expanding the above result in powers of — 1/Loge yields 



me 



ttA 2 



Expansion (vi) implies that 



Log a 



Loge 



Loge 



vi 



(— 

\Loge 



[u(x)f 



0) 



(noting that = f° r this example) which directly contradicts expan- 

sion (v), i.e. (38), except possibly in the special case a — 1. 



References. The authors of iRocha de Faria et al.l (120071 ) were apparently 
not aware of recent references directly related to their work , in particular stud- 
ies concerned with small-defect asymptotic expansions (e.g. lAmmari and Kang . 
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2004 ; IVogelius and Volkovl . I200CH ; IVolkovl . 120031 . and works cited therein) and 



with the top ological derivative for 3-D in the context of scalar and elastic wave 



propa gation (IGuzina and Chikichevi . 120071 ; iBonnetl . l2006bl ; iGuzina and Bonnet 
2006h . 
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